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We consider the initial value problem for a massless scalar field in the Schwarzschild geometry. When constructed using a 
complex-frequency approach the necessary Green's function splits into three components. We discuss all of these in some 
detail: 1) The contribution from the singularities (the quasinormal modes of the black hole) is approximated and the mode-sum 
is demonstrated to converge after a certain well defined time in the evolution. A dynamic description of the mode-excitation 
is introduced and tested. 2) It is shown how a straightforward low-frequency approximation to the integral along the branch 
cut in the black-hole Green's function leads to the anticipated power-law fall off at very late times. We also calculate higher 
order corrections to this tail and show that they provide an important complement to the leading order. 3) The high-frequency 
problem is also considered. We demonstrate that the combination of the obtained approximations for the quasinormal modes 
and the power-law tail provide a complete description of the evolution at late times. Problems that arise (in the complex- 
frequency picture) for early times are also discussed, as is the fact that many of the presented results generalize to, for example, 
Kerr black holes. 

can be approached "analytically". This can lead to an 
improved understanding of the underlying physics and 
information that can be extremely difficult to infer from 
purely numerical data. 

The problem we consider here is in many ways an old 
one. The evolution of a test-field in a black-hole back- 
ground was first considered more than 25 years ago . 
The general features of such an evolution are well-known 
. The "response" of the black hole — as seen by a dis- 
tant observer — can be divided into three components. 
Radiation emitted "directly" by the perturbation source 
will dominate at early times. This radiation depends on 
the exact character of the initial field. In contrast, the 
late-time response depends mainly on the parameters of 
the black-hole. The exponentially damped oscillations of 
the black-hole quasinormal modes carry a considerable 
part of the total radiated energy in many astrophysical 
processes (such as gravitational collapse) |5[^,0,§- Fi- 
nally, the wave field falls off with time according to a 
power law at very late times |^,|o| . 

The initial-value problem for black-hole perturbations 
have been considered by several authors. In an impressive 
study. Leaver discussed both the excitation of quasi- 
normal modes and the nature of the power-law tails. The 
quasinormal-mode problem was later considered by Sun 
and Price and also by the present author Q . quasi- 
normal modes in physics studied by Gundlach, Price and 
PuUin |l|,|l| and Ching et al. 

Even though the problem is far from new, there are 
several reasons why it need be investigated further. Al- 
though the response of a black hole to an impinging 
wavepacket will almost exclusively be dominated by the 
slowest damped quasinormal modes — and present meth- 
ods can reliably account for the excitation of these modes 
— several questions remain. For example: What is 
the role of the highly damped modes? It is known that 
an infinite number of quasinormal modes exist for each 
radiative multipole £ jl^jl^, but our understanding of 
the role of the higher overtones is rather poor. In fact. 



I. CAUCHY'S PROBLEM FOR PERTURBED 
BLACK HOLES 

This paper concerns the evolution of a test-field (be it 
scalar, electromagnetic or a perturbation of the gravita- 
tional field itself) in a spacetime that contains a black 
hole. That is, we consider the problem that is associ- 
ated with Cauchy in the framework of general relativity. 
Because of the inherent nonlinearity of Einstein's theory 
this problem is generally not amenable to analytic calcu- 
lations. But if the wave-field is sufficiently weak that its 
contribution to the spacetime curvature can be neglected 
the evolution equations reduce to a wave equation with 
a complicated effective potential. This is the realm of 
black-hole perturbation theory in which the initial- 
value problem can be approached by "standard" methods 
|,|. The purpose of the present work is to contribute 
a more detailed understanding of the many intricacies 
associated with the evolution of a weak wave-field in a 
black-hole geometry. 

One can argue that this kind of discussion is of little 
importance to physics. It may seem obvious that much 
relevant information will be lost when the equations of 
general relativity are linearised. But it turns out that the 
perturbation approach provides surprisingly accurate re- 
sults in many situations. An interesting example of this is 
the case of two colliding black holes . This obviously 
does not mean that the linear equations render a fully 
nonlinear approach useless. It would be truly surprising 
if no new phenomena were to be unveiled by detailed non- 
linear calculations, but linear studies provide important 
benchmarks against which such fully nonlinear, numeri- 
cal calculations can (and should) be tested. Also — and 
of equal importance — is the fact that the linear problem 
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it is not at all clear whether the mode-sum is conver- 
gent or not m . Our understanding of the power-law tail 
is also somewhat unsatisfactory. The leading behaviour 
has been calculated in different ways HJl^ , but the resul- 
tant formulae are only truly useful at very late times. In 
a typical evolution scenario there is a considerable time- 
window in which the signal is no longer dominated by the 
quasinormal modes, but the leading order power-law tail 
has not yet taken over. Is it possible to derive a "higher- 
order" tail expression that describes the evolution ade- 
quately for the intermediate times? These questions (and 
several others) are addressed in the present paper. 



II. THE PROBLEM AND A FORMAL SOLUTION 

A. A massless scalar field in the Schwarzschild 
geometry 

In order to make the presentation clear we have chosen 
to specialize this investigation to the case of a massless 
scalar field and Schwarzschild black holes. This is, of 
course, a model problem since no scalar fields have yet 
been observed in Nature. But this does not mean that 
our results are of restricted value. On the contrary: Be- 
cause the equations that govern other perturbing fields 
(such as an electromagnetic testfield or a gravitational 
perturbation of the metric) are similar to the one for a 
scalar field Q, the results presented here are easily ex- 
tended to all other relevant cases. Furthermore, as will 
be discussed in section VI, it seems likely that many of 
the present results can be adapted also to the case of 
rotating black holes. 

In the background geometry of a Schwarzschild black 
hole a massless scalar field evolves according to 
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Because of the underlying spherical symmetry it is mean- 
ingful to introduce the decomposition 



(2) 



where Y^m are the standard spherical harmonics. The 
function Ui{r^,t) then solves the wave equation 



where the effective potential is 
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and M is the mass of the black hole (we use geometrized 
units c = G = 1). The "tortoise" coordinate is defined 

by 



d 



= l- 



2M 



dr 



(5) 



Let us now suppose that we are given a specific scalar 
field at some time (we will use t = 0), and that we want 
to deduce the future evolution of this field. That is, we 
require a scheme for calculating (for each I) U£{r^:, t) once 
we are given u^(r»,0) and 9tit^(r»,0). This problem is 
typically approached via a Green's function. 



B. Tlie blacl4-hole Green's function 

It is well-known that the time-evolution of a wave-field 
ue{r^,t) follows immediately from 



G{r^,y, t)dtUi{y, 0)dy 
dtG{r^,y,t)ue{y,0)dy , 



(6) 



for t > (we will discuss the appropriate limits of inte- 
gration in section IIIC). The (retarded) Green's function 
is defined by 



G{r,,y,t)^S{t)S{r,~y) , (7) 



together with the condition G{r^,, y,t) — for < < 0. 

To find the Green's function G{r^,,y,t) is our main 
task. Once we know it, we can study the evolution of 
any initial field by evaluating the integrals in (^. 

The first step in finding G{r^, y, t) consists of reducing 
(0) to an ordinary differential equation. To do this we 
use the integral transform B 



G(r*,y,w) 



G(r,,y,t)e'"*dt . 



(8) 



This transform is well defined as long as Im lu > 0, and 
the corresponding inversion formula is 
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G{r,,y,t) = — G{n,y,u;)e-'^'du , 

^'^ J —oc+ic 



(9) 



where c is some positive number. 

The Green's function G{r^, y, uj) can now be expressed 
in terms of two linearly independent solutions to the ho- 
mogeneous equation 



(10) 



The two required solutions are defined by their asymp- 
totic behaviour. The first solution corresponds to purely 
ingoing waves crossing the event horizon: 



^out(t^)e'"'-- + An(w)e- 
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and the second solution behaves as a purely outgoing 
wave at spatial infinity: 



uf{r^,uj) r ^ ^^^^^ 



(12) 



Using these two solutions the Green's function can be 
written 



(13) 

Here we have used the Wronskian relation 



2iujAin{uj) . (14) 



C. Using complex frequencies 

The problem can now, in principle, be approached by 
direct numerical integration of (|lO|) for (almost) real val- 
ues of ll> and subsequent inversion of (|^). This approach 
should lead to reliable results and an accurate represen- 
tation of the evolution, as long as some obvious care is 
taken in each step. A multitude of examples of this ap- 
proach can be found in work relating to particles orbiting 
black holes (see |^] for an exhaustive list of references). 
For the evolution of a test-field, when we want to explain 
why different features seen in the emerging waves arise, 
it may be useful to follow an alternative route, however. 
An approach that often proves useful when one wants to 
isolate the behaviour of a Green's function in different 
time intervals is based on bending the integration con- 
tour in (^) into the lower half of the complex w-plane. 
This is the approach that we will follow here. 

What do we expect to learn by analytically continu- 
ing the Green's function in this way? First of all, it is 
well known that G(r, , y, uj) has an infinite number of dis- 
tinct singularities in the lower half of the w-plane. These 
correspond to the black-hole quasinormal modes and oc- 
cur at frequencies for which the Wronskian W{'w) van- 
ishes. That is, for a quasinormal mode the two solutions 
and m"'' are linearly dependent. To determine the 
quasinormal-mode frequencies is not a trivial task, but 
several accurate methods have been devised |l^j2^,^,^ . 
The mode frequencies do not, however, contain all the in- 
formation that is required to evaluate the Green's func- 
tion. While it is formally straightforward to use the 
residue theorem to determine the mode-contribution it 
is, in practice, non-trivial to evaluate the resultant ex- 
pressions. One must be able to approximate the eigen- 
function associated with each quasinormal mode. 
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FIG. 1. a) Schematic description of a black holes response 
to initial data of compact support. The directly transmitted 
wave (from a source point y) arrives at a distant observer (at 
r,) roughly at t — + y = 0. The black holes response, that is 
dominated by quasinormal mode ringing, reaches the observer 
at roughly t — r* — ?/ = 0. At very late times the signal falls 
off as an inverse power of time. This power-law tail arises 
because of multiple backscattering off the spacetime curva- 
ture, b) Integration contours in the complex frequency plane. 
The original inversion contour for the Green's function lies 
above the real frequency axis. When analytically continued 
in the complex plane this contour can be replaced by the sum 
of 1) the quasinormal modes [the singularities of G{r,,y,uj); 
the first few are represented by crosses in the figure] 2) an 
integral along the branch cut (a thick line along the negative 
imaginary uj axis in the figure), that leads to the power law 
tail, and 3) high frequency arcs (that one would expect vanish 
at most times, but they should also lead to rougly "flat space 
propagators" at early times). 

In the complex-frequency picture the late-time power- 
law tail is associated with the existence of a branch cut 
in u"^. This cut is usually placed along the negative 
imaginary w-axis. It has been demonstrated that the 
behaviour at very late times can be obtained from a low- 
frequency approximation of the integral along the branch 
cut 1^ . As regards the radiation that reaches an observer 
more or less directly from the source it has been suggested 
Q that it can be associated with the large-frequency arcs 
that are required to "close the contour" in the complex 
tj-plane (see Figure |l|) . One can argue that this should be 
the case in a handwaving way: For large frequencies the 
Green's function limits to the familiar fiat-space prop- 
agator As yet there are no detailed studies of the 
high-frequency problem, however. 



D. The asymptotic approximation 

In this paper we want to pursue the problem analyt- 
ically as far as possible. This means that we will often 
prefer a simplifying approximation over a less transpar- 
ent numerical calculation. The hope is that this will lead 
to a reasonably accurate description of the evolution, and 
at the same time provide better insight into the underly- 
ing physics. Once one has acquired this understanding it 
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will be meaningful to perform a more accurate analysis. 

In this context, a useful approximation follows if one 
assumes that spacetime is essentially flat in the region 
of both the observer and the initial data (that should 
be of compact support). We consequently assume that 
i) the observer is situated far away from the black hole. 
This means that r*/M >> 1 in (p^). ii) the initial data 
has considerable support only far away from the black 
hole. This implies that only the region where y/M >> 1 
contributes significantly to (^). 

To make life easier we will also assume that the initial 
data has no support outside the observer (only y < i\ 
are relevant). With all these restrictions the frequency- 
domain Green's function (|l^ ) simplifies to 



G{r^,y,uj) « 



A 



A, 



out ^iuj(r,+y) 



(15) 



In the following we will refer to this as the "asymptotic 
approximation" since it follows when we use the large- 
argument asymptotics for u™ and u"^ in (^3|). The use- 
fulness of this approximation should be obvious. 



III. QUASINORMAL MODES 

A. Mode-contribution to the Green's function 

As already mentioned, the quasinormal modes corre- 
spond to complex frequencies (uJn) for which the Wron- 
skian W{w) vanishes. This means that Ain{ujn) — and 
consequently it is useful to define a quantity a„ by 



Ain{uj) ^ {lO - UJn)an , 



(16) 



in the vicinity of the mode. Then it follows from the 
residue-theorem (and the fact that modes in the third 
and fourth quadrant are in one-to-one correspondence, 
see Figure 4.4 in [^) that the total contribution from 
the modes to the time-domain Green's function can be 
written [H 



C?'3(n,y,i) =Re 



Here we have defined 



Bn — 



n=0 



-iijji (t- 



^out(^n) 



(17) 



(18) 



We have also used the asymptotic approximation (15), 
and the sum is over all quasinormal modes in the fourth 
quadrant. That this expression provides an accurate 
representation of the mode-excitation has already been 
demonstrated ||^,|l3|,^ . Typical results obtained using 
( |l7| ) are shown in Figures 3a-c in Q]. 



B. Convergence of the quasinormal-mode sum 

When one evolves a test-field in the Schwarzschild ge- 
ometry one typically finds that the response of the black 
hole that is associated with the data at y [cf. (|^)] reaches 
the observer roughly when t — r* — ?/ = 0. This is not too 
surprising: The slowest damped quasinormal modes can 
be associated with the peak of the effective potential . 
Hence, one would expect the response to follow once the 
specific part of the initial data has had time to reach the 
peak of the potential (roughly at — 0) and then travel 
back to the observer. 

An obvious question concerns the convergence of the 
quasinormal-mode sum (p^. At what times (if any) will 
the sum be convergent? Previous evidence for gravita- 
tional perturbations and the first seven modes |^ sug- 
gests that the mode sum is convergent at late times, but 
fails to provide a lower limit of t at which this conver- 
gence starts. As far a late times are concerned it is easy 
to convince oneself that the sum should converge: Two 
consecutive terms in ( |l7|) yields the ratio 



B. 



n+l 



Br, 



-i(w„+i-cj„)(t-r*-y) 



Now we know that 0,|8| 



n+ 



iM « w„M - i/A 



oo 



(19) 



(20) 



Assuming that the term in the square brackets remains 
of order unity (say) it follows that the magnitude of the 
ratio of successive terms in the mode sum behaves asymp- 
totically as exp[— (t — r^, — y)/4]. This implies that the 
sum will surely converge for t — — y ^ 0. 

But this argument relies on the terms in the square 
brackets of (^9|) behaving in a certain way. Does it hold 
in practice? To test this we have used the approximate 
formulae derived in Q to obtain _B„ for the first 200 
modes of a scalar field (and £ ~ 0,1 and 2). These re- 
sults are illustrated in Figure |^. The data here is not ex- 
pected to be very accurate for the highly damped modes. 
Nevertheless, the trend is clear. The magnitude of Bn 
decreases monotonically for large values of n. Moreover, 
one finds that successive terms have opposite signs. This 
suggests a stronger convergence than the expected one: 
The mode sum will converge also at t — — ?; = . 
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FIG. 2. The absolute value of the terms in the mode sum 
&t t — r, ~ y — for £ — 0,1 and 2 are shown as a function 
of the mode index n for the first 200 modes. The data is 
obtained using approximate phase-integral expressions. 

C. Dynamic mode-excitation 

Previous investigations of this problem |^,|l^ were to 
a certain extent marred by what can be called the "tim- 
ing problem" . Assume that the initial data consists of 
a "mountain" close to the black hole and a tiny "pim- 
ple" far away. It was then found that the "pimple" leads 
to a much larger mode-excitation than does the "moun- 
tain". This is, of course, contrary to our expectations. 
Fortunately, it is also wrong. 

The "timing problem" arises when one tries to asso- 
ciate a given set of initial data with a constant "excita- 
tion strength" of each quasinormal mode. This would be 
a useful approach for many familiar oscillating systems, 
such as a vibrating string. But this approach is probably 
only meaningful when the modes of the system form a 
complete set. In the black- hole case it is well known that 
the quasinormal modes are not complete (one must also 
take account of the branch cut integral). As we shall see, 
it makes more sense to consider the quasinormal-mode 
excitation as a dynamic process. 

For example, one would not expect the quasinormal 
modes to be (considerably) excited until after the rele- 
vant feature in the inital data has scattered off the poten- 
tial barrier that surrounds the hole. In our previous ex- 
ample this means that the mode-excitation that arises be- 
cause of the "mountain" in the data will be relevant ear- 
lier than that associated with the "pimple" . Roughly, the 
modes should be excited when the relevant data reaches 
the peak of the effective potential (r* r:! 0). 

It is not very difficult to obtain a dynamic descrip- 
tion of the mode-excitation. In fact, we need only ensure 



that the Green's function respects causality. The basic 
"mistake" of the previous studies |§,|3| was to perform 
the integral over the product of the initial data and the 
Green's function G{r,,,y,uj) before inverting the integral 
transform. If one instead uses (0) the confusion of the 
"timing problem" can be avoided. 

To ensure causality we must ensure that the Green's 
function vanishes when t — + y > 0. An observer (at 
r*) should simply not see anything until after a signal 
travelling at the speed of light has been able to reach 
him/her from the relevant source point (y). In the evo- 
lution equation (|^) this translates into a lower limit of 
integration y = — t. 

For the quasinormal modes, one can also argue for the 
use of an upper limit of integration: If we use (|5|) it is 
clear that the quasinormal-mode part of the signal arises 
from 

G{r.,y,t) --^ I i^e— . (21) 

Intuitively, one would not expect it to be meaningful to 
close the integration contour C in the lower half-plane 
unless t — r^ — y > Q. At earlier times it seems unlikely 
that the contribution from the necessary high-frequency 
arcs will vanish. We will consider this issue in more detail 
in section V. For now we are content to deduce that this 
introduces an upper limit of integration y — t — r^ in the 
evolution equation (^. 

As a simple example of the implications of this discus- 
sion we consider the static initial data 

u(r,, 0) = exp[-0.05(r,/M - 400)^] , (22) 
9tu(r,,0)=0. (23) 

It is straightforward to multiply this data with 
G'^{rt:,y,t) from ( |l7| ) and then integrate (numerically) 
from y = r<f — t to y = t — r^. The result of this calcula- 
tion is displayed in Figure ^. Here we have assumed that 
the observer sits at — 500 Af. This means that, accord- 
ing to the "dynamical" description the modes should not 
be present in the signal before t k, 950Af . 

From the data presented in Figure we can conclude 
that our approximation for the quasinormal mode ringing 
is quite accurate. It is conceivable that the remaining dis- 
crepancies at later times would disappear if we used the 
true mode-functions instead of the "asymptotic approxi- 
mation" . The most interesting part of Figure ^ concerns 
the early times. It is notable how nicely the idea of a "dy- 
namic" mode excitation works in practice. This allows us 
to discuss the relevance of the high-order modes in more 
detail than what has been possible before |^,^,^. Since 
they are rapidly damped, one would expect the highly 
damped modes to be relevant only at early times. As 
can be seen in Figure ^ this is, indeed, the case. While 
the slowest damped mode represents the signal well af- 
ter (say) t ~ 965M, a much better approximation (valid 
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from (say) t « 953M) is obtained by using the sum of the 
first six modes. Should we include further modes in the 
sum this trend continues, but it becomes hard to distin- 
guish any improvement. This evidence clearly supports 
the notion that the mode-sum converges for all times af- 
ter t — — y = 0. 



0.15 




980 
time 



1020 



FIG. 3. Comparing the field obtained through evolving the 
scalar wave equation to the approximate contribution from 
the quasinormal modes. The graph shows data for £ = 2. 
The solid line represents the true scalar wave, while the two 
dashed lines are for a) the slowest damped quasinormal mode 
and b) the sum of the first six modes. 



IV. THE LATE-TIME POWER-LAW TAIL 

It is by now well-known that the quasinormal-mode 
ringing is followed by a power-law tail at late times. This 
feature was first found in the seminal work of Price |^,|lO| . 
Physically, the tail arises because of backscattering off 
the slightly curved spacetime in the region far away from 
the black hole [Q. This means that the tail will not 
depend on the exact nature of the central object. Thus, 
a neutron star of a certain mass will give rise to the same 
tail as would a black hole of the same mass. 

Mathematically, it has been demonstrated that the 
power-law tail can be associated with the branch cut in 
the complex- frequency Green's function Using this 
fact, the exact form of the leading order tail has been 
calculated in different ways. But some questions still 
remain. The most important one concerns the black- 
hole response at intermediate times — after the quasi- 
normal modes have died away, but before the leading 
tail term accurately represents the evolution. Is it pos- 
sible to extend existant calculations in such a way that 
one can approximate the evolution also at these inter- 
mediate times? Previous calculations are also somewhat 



involved; it would be nice to have a simpler and more 
direct calculation of the tail effect for black holes. 

One can argue that the late-time behaviour should fol- 
low from the low-frequency contribution to the Green's 
function. Basically, the effective black-hole potential in 
( pO| ) will be so small for large values of r that only low- 
frequency waves will be affected by it. Hence, a low- 
frequency approximation to the black-hole equation will 
be useful. In this section we obtain such an approxima- 
tion, and use it to study the detailed behaviour of the 
power-law tail. 



A. A low-frequency approximation 

Let us begin by introducing a new dependent variable 

-1/2 

Then equation (|10|) for the scalar field becomes 



(24) 



1 



2M 

r 



1 



2M\ £{£+!) 
r J 



— 

(25) 



dr'^ 



t/- w 



(26) 



rtp = 2M 



AM 



2, ,,2 



-2-1 



(27) 



(m) ^'^'^ 



(28) 



z— + {2l + 2- z)— - (^ + 1 - 2iujM) 
dz^ dz 



(t> = Q . (29) 



= 



As already mentioned above, it is easy to argue that we 
only need a large r approximation to account for the low- 
frequency response. Thus, we expand (|2^) as a power 
series in r/M. This leads to 



The question is whether this is a useful approximation 
psf . That is, is the assumption that we only need large 
r justified? Fortunately, this is a simple thing to test. 
From ( p^ ) we find that the outer turning point of classical 
motion is located at 



At rtp one would expect a wave with frequency lj to be 
scattered by the effective potential. Since rtp — s- -l-oo as 
ljM ^ it is clear that the approximate equation ( |2^ ) 
can be used with confidence for low frequencies. 
Let us now introduce |29| 



where z = —2iujr. Then, it follows that (j) should be a 
solution to the confiuent hypergeometric equation 
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It also follows that two basic solutions to the black-hole 
problem can be written (remember that r >> 2M) 

iif ^ A [j;^y^^ e*"''M(^ + 1 - 2iujM, 21 + 2, -2itjr) , 

(30) 

and 

-up _ 



M 



u^'P = S — e*"'^C/(£ + 1 - 2iujM, 2£ + 2, -2iujr) 



(31) 



where A and B are normalisation constants. The func- 
tions M{a, b, z) and U (a, 6, z) represent the two standard 
solutions to the confluent hypergeometric equation [2^ . 

One nice feature of this approximation is that it is 
obvious that there will be a cut in u^^. We simply use 
the standard result |2^ that, for n an integer, 

(-!)"+! 

U{a,n+ 1, z) = -Mia, rt + 1, z) Inz + 

n\V (a — n) 

+ single-valued terms . (32) 

From this it follows immediately that 

U{a, n + 1, ze^") = U{a, n + 1, z) 

+ 2m^—^ -M{a, n+l,z) . (33) 

n\l (a — n) 

We will now use this result to evaluate the effect of the 
branch cut in the Green's function. 



B. The late-time tail 

We have seen that there will necessarily be a branch cut 
in the black-hole Green's function G{r^,y,ijj) (or more 
specifically, in the solution u^')- To arrive at the detailed 
power-law tail we need to consider the effect of this cut. 
It is easy to see that the contribution from the branch 
cut follows from the integral [cf. (^)] 



G^(r-*,2;,i) = 

1 p—iOQ 

u"''(r*, w) 



(34) 



W{uj) 

After some straightforward steps this can be rewritten as 
G^(^*,y,0 = 

2^/ 



M^[a;P(r„L^),u;P(n,c^e2")] 



Using the low-frequency approximation that was de- 
scribed in the previous section we can show that 

and since there is no cut in u^" it is easy to see from ( |3^ ) 
that 

Wicue^"") = W{lo) . (37) 
To build the Green's function we also need 
I¥[w;P(r„c^),u;P(r„c^e2")] = 



B 



, (-l)-^-i27r(2^)-2^-i 



(38) 



T{-e - 2iujM)r{e + 1 - 2iLuM) 

If we use these results together with the approximation 

uj," « A [■^Y'^^ e''^''M{(. -f 1, 2^ 4- 2, -2iujr) = 

^ A{2i+l)\\{ojM)-' n{ur) , (39) 

we get 

n—ioo 

JQ 

(40) 

To obtain this result we have assumed that the asymp- 
totic approximation is valid. Specifically, we have re- 
placed r by This is not at all a necessary step, but 
it simplifies the comparison to our numerical evolutions 
of ^ which were carried out using r, as independent 
variable. Moreover, if r is sufficiently large we will only 
introduce a small error. In the specific case of ( p3| ) the 
error introduced in this way will be smaller than 2.5 % . 

As long as we are only interested in the leading order 
behaviour at very late times, we can assume that cuy < 
Lur^ « 1 Thus, using the standard power series 
expansion for je{z) we arrive at the final formula 



G^(r„y,t) = (-1)^ 



_i (2^ + 2)! AM{r,yY+^ 
[(2^+l)!!]2 t2^+3 



(41) 



(35) 



This result is identical to that obtained by Leaver 
This is not surprising since our derivation is just a sim- 
plification of Leaver's approach. Our result also agrees 
with that of Ching et al [ p^ , but in that case one must do 
some additional calculations to obtain the explicit result 
for the black-hole case. 

We have not yet contributed much new information. 
The greatest merit of the above derivation is its simplic- 
ity. The origin of the branch cut is clear and its contri- 
bution to the Green's function follows painlessly. This 
is in remarkable contrast to, for example, the work of 
Ching et al p6|], where the result follows after a truly 
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involved analysis. On the other hand, the formulae ob- 
tained by Ching et al. are valid for a large class of poten- 
tials. Hence, their work shows that the tail-phenomenon 
is a generic feature of many problems in wave scattering. 

Simplicity is not the sole advantage of the present ap- 
proach, however. It turns out that we can easily do the 
full integral in (^). As long as t — r, — y > we can em- 
ploy equation (6.626) from Gradshteyn and Ryzhik [ p9[ . 
This leads to a higher-order result 



M 



i-lY 



'{2£ + 2m + 2)l 



oo m 

EE 

nUm - n)U2£ + 2n+ 1)!!(2£ + 2m-2n+l)\\ 



^e+2m~2n+l e+l+2n 
V — - 

which is remarkably simple. 



(42) 



also clear from Figure ^ that if one includes several terms 
in ( ^ ) one arrives at an approximation that takes over 
from the quasinormal-mode ringing in a natural way. 

Figure ^ shows that there will be a considerable time 
window in which the contributions from the quasinormal- 
modes and the higher order tail are of the same order of 
magnitude. It seems reasonable to assume that this may 
lead to interference effects that can be distinguished in 
the evolution of the field. That such effects are present 
is clear from Figure ||. It is also clear that a combination 
of the quasinormal modes and the higher order tail-sum 
provides a good representation of the evolution through- 
out the transition between the typical times when either 
term dominates the signal. Thus the approximations dis- 
cussed so far can be combined to estimate the signal for 
all times t — — y > 0. 
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FIG. 4. Comparing the field obtained through evolving the 
scalar wave equation to the approximate contribution from 
the power-law tail. The graph shows data for £ — 2 using 
a logarithmic scale. The solid line represents the true scalar 
wave, while the three dashed lines are for a) the leading order 
power-law tail b) the tail approximation including the first 
two terms and c) using the first 11 terms in the tail-sum. 

We have thus managed to extend previous work to in- 
clude higher order corrections to the power-law tail. The 
question is to what extent such a result is useful. Inter- 
estingly, it turns out that the higher order terms play an 
important role. Essentially they allow us to extend the 
validity of the tail-approximation to much earlier times. 
A typical result — obtained for the initial data given in 
( p3| ) — is shown in Figure ^ Although it is not clear 
from Figure ^, it can be verified that the lowest order 
tail-term is a reasonable approximation to the true evo- 
lution at very late times. But the improvement achieved 
by including the first two terms in (E4) is impressive. It is 
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FIG. 5. Comparing the field obtained through evolving the 
scalar wave equation to the approximate contribution from 
the slowest damped quasinormal mode and the first 11 terms 
in the sum for the power-law tail. The graph shows data for 
£ — 2 using a logarithmic scale. The solid line represents 
the true scalar wave, while the dashed line represents the ap- 
proximation. It is easy to distinguish eflFects of interference 
between the mode and the tail-terms (whenever the ringing 
diflFers from pure expontential damping at a constant oscilla- 
tion frequency). 



V. HIGH FREQUENCIES 

In the complex-frequency approach it is typically as- 
sumed that the contribution from the required arcs at 
|cl;| ~* 00 is irrelevant. That way, the original contour- 
integral can often be replaced by a mode-sum such as 
(|l7|). In the black-hole problem the situation is, of course, 
complicated by the branch cut in the Green's function. 
But as we have seen, the contribution from this cut can 
readily be approximated (at least at late times). To com- 
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plete the study of the initial value problem we now focus 
our attention on the high-frequency problem. 

Although it is easy to give a handwaving argument 
that suggests that the large-frequency arcs should not 
contribute significantly to the Green's function (when 
\ojM\ is very large one would not expect the details of 
the effective potential to matter much; a high-frequency 
wave will propagate almost as in flat space) the issue has 
not been studied in much detail previously. A better un- 
derstanding of the high-frequency problem is useful for 
several reasons: i) We obviously should confirm that our 
expectations of a vanishing contribution to the Green's 
function holds, ii) The high frequencies may hold the 
key to the response of the black hole at early times 
Hence, it is plausible that an understanding of the be- 
haviour for high frequencies will yield a handle on the 
inital part of the signal that reaches an observer. 



A. An approximation for high frequencies 

An approximation that is relevant for high frequencies 
can be obtained in the following way : When \ujM\ be- 
comes large the equation that governs the scalar field (|2£^ 
limits to the confluent hypergeometric equation |28[ 



az^ az 



= . (43) 



Here we have used 
^ \2M 



ij^exp [iuj{r -2M)](j) 



z = -2iuj{r-2M) 



and ^ is given by 



2iujM 



(44) 



(45) 



(46) 



As in the case of low frequencies it is easy to use stan- 



dard formulae from |27| and show that the two basic so- 
lutions that we require to build the black-hole Green's 
fimction can be written 



xM(l/2 - 4iwM, 1 - AiujM, -2iuj{r - 2M))e'"('^"'^*^) 



-11 X 



and 



-.up _ 



xJ7(l/2 - AiujM, 1 - AiujM, -2iu{r - 2M))e" 



-4iLuM 
1) X 



_ r(l - 4itjM)(4itjM)-l/2+4»-'Mg-4».^M 

^out — 1= , (49) 



and 



r(l - 4iwAf)(4iwM)-i/2e*'^/2 
T{l/2- iiujM) 



(50) 



After using large argument approximations for the T 
functions we get 



and 



(51) 



(52) 



Hence, a high-frequency approximation to the "reflection 
coefficient" of the black hole is 



n 



A. 



(53) 



This result agrees with our expectations (from for exam- 
ple the WKB method): For very large frequencies the 
reflection caused by the black-hole potential barrier will 
be exponentially small. 

We can also use this approximation to approximate 
the very high overtones of the black hole. Recall that 
the quasinormal modes follow from A-^ = 0. Then, it is 
trivial to use (^0|) and show that modes should be located 
at M 



(54) 



This approximation yields the correct damping rate for 
the high modes, but it fails to reproduce the small con- 



stant real part that each mode should have |17 1 



B. High-frequency Green's function 

We now want to use these results to discuss the possi- 
ble contribution to the black hole Green's function from 
the required arcs at \ujM\ 00. To do this we assume 
that the asymptotic approximation from section HD is 
appropriate (we will discuss alternatives to this later). 
Then we clearly get (in the right half-plane) 



(47) 



G°°(r„2/,c.)«--^ 
2iuj 



(55) 



(48) 



From the asymptotic behaviour of the confluent hyper- 
geometric functions ||2^ it follows that (in the right half 
of the complex a;-plane) 



For obvious reasons it makes sense to study the two terms 
in the square brackets separately. 

For each term we require an integral of form [ cf. (0)] 



(56) 
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where C is a large-frequency quarter circle in either the 
upper or the lower right half of the w-plane (the calcula- 
tion here must be complemented by similar formulae for 
the left half- plane). Now, as long as 



, as 



oo 



(57) 



(which is certainly true here) the integral / vanishes i) 
in the upper half-plane for r < 0, and ii) in the lower 
half-plane for r > (Jordan's lemma). 

For the specific case of the black-hole Green's function 
this means that 



G{r,,y,t) = 0, for 



t-r^+y <0 
t — — y > 



(58) 



This is a nice result because it shows that a combination 
of the quasinormal modes and the contribution from the 
branch cut in G(r, , y, uj) should form a complete descrip- 
tion of the wave evolution after t — r^, — y = 0. This 
confirms the result of the previous section. It also shows 
that the evolution should be causal: No signal will reach 
the observer before t — r,, -\- y = Q. But the situation for 
the intermediate times is not clear. 



C. Approximating the signal at earlier times 

We have seen (cf. Figure ||) that the combination of 
the quasinormal modes and the higher order tail sum pro- 
vides a complete description of the evolution of an initial 
wave field after a certain time. The only remaining ques- 
tion is whether we can approximate the evolution ade- 
quately also at earlier times. Unfortunately, this turns 
out to be harder than what one might expect. 

To discuss this issue we consider the high-frequency 
Green's function (|55|). For intermediate times ?'* — ?/< 
t < -\- y one would naively want to use a contour in 
the lower half-plane for the first term in (^^ while at 
the same time using a contour in the upper half-plane 
for the second term. But although consistent with our 
discussion of quasinormal modes in section 111 this is not 
a useful approach. The main reason is that we should 
not treat the two terms in ( |5^ ) independently. This does 
not mean that our previous discussion is flawed: Our 
results still hold because after t — r* — ?/ > one would 
certainly want to close both integration contours in the 
lower half-plane. But problems arise when we want to 
consider earlier times. 

Before discussing this problem in more detail it may 
be helpful to illustrate the features in the evolution that 
we want to describe. In Figure ^ we show the initial 
signal that reaches the observer (at t — + j/ = 0) in 
the case of (|2^). It is clear that this signal is more or 
less the direct transmission that one would expect in flat 
space. But there is one important difference that is hard 



to distinguish in Figure ^: After the initial pulse follows 
a tiny wake. A description of the early time behaviour 
must yield this effect. Another part of the signal that 
one would like to describe is the "reflected" pulse that 
reaches the observer slightly before t — — y = 0, i.e., 
before the onset of quasinormal ringing in Figure ^ 

So why is this problem difficult? Basically, there are 
two possible routes and both forces us to deal with for- 
mally singular terms. The first possibility is to close the 
integration contour in the upper half-plane at all times 
before t — — y = 0. If we do so it is clear that the 
first term in (|5^) will be singular, but all other contri- 
butions to the Green's function will vanish. The second 
option is to close the contour in the lower half-plane as 
early as t — + y = 0. Then we find that there will be 
three divergent terms (that must balance each other in 
some magic way): We know that the quasinormal-mode 
sum will diverge, and the same is true also for the inte- 
gral along the branch cut and the high-frequency arc [the 
second term in (|55|)]. Clearly, this second alternative is 
the less attractive one and we will take a few steps down 
the first route here. 
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FIG. 6. Comparing the field (for = 2) obtained through 
evolving the scalar wave equation to an approximation that 
corresponds to propagation in flat space. The initial signal 
that reaches the observer is shown. After the, essentially un- 
changed, initial shape follows a tiny wake. This wake will be 
more pronounced for data that has support closer to the black 
hole. The solid line represents the true scalar wave, while the 
dashed line represents the approximation. 

Let us consider only the leading order. As already 
pointed out, the only contribution to the Green's function 
that does not vanish for Vi^— y<t<r^-\-y comes from 
the first term in (|55|). Thus, we need to evaluate 



G{r,,y,t) 



-I- 

47ri Jc 



-iuj{t — r^ +y) 



-dtj 



(59) 
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where the integration contour C is a semi-circle in the 
upper half-plane. Now the integrand can be identified as 
the Laplace-transform of the step-function. Specifically, 
we get 

G{r,,y,t)^^H{t-n-y) . (60) 

This is, of course, exactly what we would get in flat space. 
To understand the added subtleties of the black-hole case 
we must pursue the calculation to higher orders. 

In principle, such a higher order calculation seems pos- 
sible. One could, for example, try to express the high- 
frequency approximations of the solutions ii^^ and uj," 
from section VA as power series in w"". That should 
yield higher order corrections to (|60|). Unfortunately, it 
seems as if one would have to keep a large number of 
terms (all?) in the resultant expression to get a useful 
answer. There may also be less obvious complications. 
Hence, we will have to return to this problem in the fu- 
ture. 

VI. EXTENDING THE PRESENT WORK 

To conclude this paper it is meaningful to discuss 
how the present work — for a massless scalar field in 
the Schwarzschild geometry — can be adapted to other, 
physically more interesting, cases. 

A. Other perturbing fields 

The simplest extension of the present work regards 
other perturbing fields. It is, in fact, trivial to show 
that all our results carry over also to electromagnetic 
and gravitational waves in the Schwarzschild background. 
Both the discussion of the tail effect in section IV and 
the high-frequency discussion in section V are valid also 
for these other test-fields. As regards the quasinormal 
modes, the characteristic frequencies (and the coefficients 
Bn in (p^)) will be different for other fields, but the ap- 
proximate phase-integral expressions |^ that we used to 
evaluate ( [l7|) can be used also for electromagnetic and 
gravitational perturbations. 

B. Without the asymptotic approximation 

Throughout this paper we used the asymptotic approx- 
imation from section IID to simplify the calculations. 
This clearly restricts the initial data in an unneccessary 
way. Fortunately, it is not (formally) difficult to general- 
ize our results in such a way that they hold also for more 
general data. First we note that the asymptotic approxi- 
mation was never used (apart from in the replacement of 
r by r* in (HQ)) in the derivation of the tail-expressions. 



Consequently, these results remain unchanged also for 
general initial data. For the quasinormal modes, we must 
use approximations of the corresponding eigenfunctions 
that remain valid for all r. That this can be done has 
already been demonstrated by Leaver [Q. The dynamic 
mode-excitation from section IIIC becomes more difficult 
to introduce in the general case. With the asymptotic 
approximation it was easy to find a time before which it 
would be meaningless to use integration contours in the 
lower half of the w-plane. In the general case, this specific 
time might not be so easy to define. This issue is inti- 
mately related to the high-frequency problem. It is clear 
that the asymptotic approximation is crucial for the dis- 
cussion in section VB. In the general case one should use 
large-w asymptotics for (^7|) and ( ^ ) that do not at the 
same time assume a large value of r. Such asymptotic ex- 
pressions are not standard but they should be possible to 
derive Q . These new approximations should indicate at 
what time it is sensible to use integration contours in the 
lower half of the w-plane, and thus imply an initial time 
for the excitation of quasinormal modes. Consequantly, 
it seems likely that all the ideas presented in this paper 
can be useful also in the general case. 

C. Rotating black holes 

We finally turn to the problem for rotating black holes. 
When the black hole is rotating both the angular and the 
radial functions that are required to describe a pertur- 
bation are frequency-dependent. In a numerical evolu- 
tion one would therefore not decouple the corresponding 
equations, and thus have to deal with a two-dimensional 
problem. The analysis of the Kerr problem is, in gen- 
eral, far more complicated than the present one. But in 
principle one would expect all the ideas discussed in this 
paper to be useful also for rotating holes. 

For example, the construction of the black hole Green's 
fmiction should be analogous to that in section II. Of 
course, the angular dependence would have to be in- 
cluded in equations like (^ . Then the quasinormal modes 
are defined exactly as in the Schwarzschild case. The 
main difference is that each Schwarzschild mode splits 
into 2£ + 1 distinct ones (for different values of to) be- 
cause of the rotation of the black hole [p^ . 

As regards the late-time power-law tail, the generaliza- 
tion to Kerr also seems straightforward. Since the issue of 
tails in the geometry of a rotating black hole has recently 
been analyzed through numerical evolutions by Krivan, 
Laguna and Papadopoulos it may be worthwhile to 
discuss this in somewhat more detail here. 

As was first shown by Teukolsky, the equations that 
govern a small perturbation of a rotating black hole can 
be written js^ 
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dr \ dr 



- 2is{r - M)K 
A 



+ AisLor — A 



i? = 0, (61) 



and 
1 d 



sindde V'"^ d^? 



a^Lu'^ cos^ I 



m 



2ms cos t 



sm 



cot^ 9 + E-s^ 



5 = 0. 



Here we have introduced 

K = {r^ + a'^)uj - am , 

and 

\ = E — s{s + 1) + a^Cij^ — 2amuj 



2aijjs cos 6 
(62) 

(63) 
(64) 



Here a < Af is the rotation parameter of the black hole, 
and s is the spin- weight of the perturbing field. The solu- 
tions to the angular equation ( |6^ ) are generally referred 
to as "spin- weighted speroidal harmonics" . 

Let us now adopt the approach of section IV and ap- 
proach the Kerr problem for low frequencies. It is suffi- 
cient to consider s > (the results for s = —1,-2 can 
be deduced via the Teukolsky-Starobinsky identities 0), 
and we can also use the fact that [p3| 



E^eie + i) + s^ + J2fniau;r 



Then we expand (pi)) for large r and find 



(65) 



dr"^ 



UIu)^ + 2isuj 



2amuj — 2isMuj 



where the new dependent variable "0 is defined by i? = 
A-('<+i)/2^^ From this equation two things follow im- 
mediately: i) For low frequencies we can neglect all fre- 
quency dependent terms in the factor of 1/r^. ii) For 
scalar waves (s = 0) we reproduce the Schwarzschild 
equation (p^). Hence, the leading-order tail calculation 
in section IV must remain valid also for scalar waves in 
the Kerr background. This confirms the main result of 
Krivan, Laguna and Papadopoulos pi| . But what about 
fields with non-zero s? Well, from p^ ) and (|6^) follows 
that we can neglect all terms that depend on the rota- 
tion of the black hole. This shows that it is sufficient 
to analyze the tail effect in the Schwarzschild limit. The 
fact that the radial equation ( |66| ) does not immediately 
agree with ( |2^ ) need not worry us too much. As Chan- 
drasekhar has shown [Q, the a ^ limit of (^l]) — 
the Bardeen-Press equation — can easily be trans- 
formed into the Regge- Wheeler equation (jl^). We thus 



have clear evidence that the leading order tail result from 
section IV will hold also for Kerr black holes. 

To derive higher order tail-corrections for Kerr de- 
mands a more involved analysis. The issue is complicated 
by the frequency-dependent coupling between ( |6l| ) and 
(|6^). Similar difficulties make the high-frequency prob- 
lem less transparent. Hence, we will not discuss these 
problems further here. What is clear is that the Kerr 
problem poses an interesting challenge, and we hope to 
be able to discuss it further in the near future. 
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